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DRAG FORCES VELOCITY GRADIENT FLOW* 


SYNOPSIS 


exploratory investigation was made the drag coefficient for circular 
cylinder influenced velocity gradient along its axis. Someof the three- 
dimensional aspects the flow mechanics were described qualitatively. Ex- 
periments indicated that the local drag coefficient varied significantly along the 
length the cylinder. For velocity gradients inthe range engineering signi- 
ficance, local drag coefficients may affected much 


INTRODUCTION 


The evaluation the hydrodynamic forces solid body flowing fluid 
isof primary importance many engineering problems. afew simple cases, 
these forces can evaluated analytical methods alone. most practical 
problems, isnecessary resort analysis and experiment. 
experimentally determined drag coefficient commonly estimating 
drag forces. 

Experimental determinations drag various shapes usually 
have been made stream with nearly constant velocity. However, many 
problems engineering importance arise which the velocity the fluid 


Note.—Discussion open until December 1960. extend the closing date one 
month, written request must filed with the Executive Secretary, ASCE. This paper 
part the copyrighted Journal the Hydraulics Division, Proceedings the Ameri- 
can Society Civil Engineers, Vol. 86, No. July, 1960. 
Presented meeting the Hydraulics Division, ASCE, Los Angeles, Calif., 

February 1959. 

Formerly, Instr. Civ. Engrg., Univ. Science-Faculty Fellow, Div. 
Hydraulic and Sanitary Engrg., Univ. California. 
Prof. Civ. Engrg., Univ. Texas, Austin, Texas. 
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stream not constant. For example, velocity variations will present the 
fluid movement the result wave motion, the fluid motion retarded 
the resistance along boundary. More specific problems, which veloc- 
ity gradient may significant, piers other structures, 
wind forces buildings stacks, resistance due piers streamsand 
channels, forces sediment particles, and surfaces irregular- 
ities boundary layer. 

Comparatively little effort has the study drag forces ina 
velocity gradient. Hoerner, (1)3 has treated two aspects the problem. 
summarizes the work Tillman, Wiegard, and others, the drag pro- 
tuberances ina turbulent boundary layer. For these surface irregularities, 
whose height less than the boundary layer thickness widthis large 
relation their height, the indications are that satisfactory estimate 
the can made assuming the drag coefficient tobe the same for 
flow with velocity gradient. effective drag coefficient based the ve- 
locity the top the protuberance can computed simply summing the 
forces successive layers accounting only for the variation the velocity. 
That is, each layer, the drag coefficient considered the same for 
two-dimensional flow. Hoerner also presents approximate data the velocity 
profiles near the ground for wind and drag coefficients for several building 
shapes. 

number papers are available, dealing with the motion solid bodies 
rotating fluids. Squire’s (2) has dealt with several aspects this 
problem and refersto group papers onthis subject presentedin the period 
from 1915-25. 

Hawthorne (6) made theoretical and experimental study the flow 
about struts and airfoils aflow with spanwise velocity gradient. His anal- 
ysis vortex development special interest, but points out not 
applicable blunt-nosed bodies. assumes that the induced spanwise velocity 
components are small and does not account for any effect the For 
these reasons his analysis, although basedon similar concepts the mechanics 
the flow, not applicable the present investigation. 


FORMULATION PROBLEM 


The variables normally consideredin studies canbe grouped into 
dimensionless equation the form 


which the drag coefficient, the drag force, thecross-sectional 
area the object taken perpendicular the flow, the mass density the 
fluid, the velocity the approaching fluid, and the Reynolds Number. 
This equation represents the average drag coefficient function the 
Reynolds Number and the shape the object under consideration. The 
application this equation two-dimensional flow problems normally assumes 
uniform velocity distribution approaching the object question. 


Numerals parentheses, thus (1), refer corresponding items the Biblio- 
graphy. 
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This investigation involved the introduction additional variables Eq. 
permit evaluation the drag forces circular cylinder with its axis 
normal the approaching flow affected linear velocity gradient the 
direction its axis. 

dimensional analysis enables gradient parameter formu- 
lated into the drag coefficient equation 


The geometrical aspects are represented the form and the cyl- 
inder diameter cylinder length, theform restricted circular 
cylinder and the ratio, held essentially constant, these two terms will 
not affect the drag coefficient. 


variables inwhich AV/AY isthe slope the velocity profile and the ratio 


the diameter the cylinder the velocity measured the midpoint the 
cylinder. the velocity distributionis uniform, the slope, AV/AY zero, thus 


causing vanish leaving the drag coefficient function the Reynolds 


Number only. 

The Reynolds Number may formed from the velocity the 
determined from the velocity profile measurements, the diameter the cyl- 
inder,and viscosity thefluid. For rangeof Reynolds numbers 
between 100,000, there practically the average drag 
coefficient with the Reynolds Number. this range Reynolds numbers, the 
shear forces the cylinder are very small relationto the pressure forces, 
and the total drag can considered being primarily form drag. 

Under the limitations previously form, D/Yo, and may 
eliminated from the functional relationship Eq. leaving the drag coeffi- 
cient function the velocity gradient parameter only. 

aid analyzing the mechanics this problem, local drag coeffi- 
cient was formulated similar the average drag coefficient, that is, 


which the local drag force, the projected area element 
the cylinder, and the local velocity the fluid approaching the element. 
The local force, may determined from integration the pressure dis- 
tribution around the cylindrical element, 


cylinder. 


Combination this local force with the corresponding measured local ve- 
locity will give the local drag coefficient. This computation, several 


which the pressure intensity any the cylinder radius,and the 
angular the cylinder measured from the upstream element the 
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elements along the longitudinal axis the cylinder, will give the variation 
the local drag coefficient along the cylinder. 

consideration the mechanics flow around the cylinder will indicate 
that the velocity gradient will components along the axis the cyl- 
inder. Along the upstream element the cylinder (the stagnation line), the 
stagnation pressure will greater the endof the cylinder, where the veloc- 
ity high, than the end velocity low. This will produce pres- 
sure gradient along the axis the cylinder andinduce flow along the stagna- 
tion line toward the low-velocity end the cylinder. similar way, the re- 
duced pressure the wake will affected the local approach velocity. 
the high-velocity endof the pressure the wake will beless than 
the low-velocity end, thus, inducing longitudinal flow toward the high-veloc- 
ity end the cylinder, the wake zone. 

The effect these longitudinal flows the local drag coefficient various 
locations along the cylinder may anticipated asfollows. the high-velocity 
end the cylinder the longitudinal flow along the stagnation line would di- 
verting fluid away and, thus, reducing the pressure intensity the upstream 
side the cylinder. the same time, the longitudinal flow the wake zone 
would forcing fluid toward the high-velocity endand thereby raising the pres- 
sure intensities the downstream side the cylinder. Both these effects 
wouldact onthe cylinder and, reduce the local 
drag coefficient the high-velocity endof the cylinder. might anticipated 
that the opposite effect would true the low-velocity end the cylinder, 
and the local drag coefficient would high. However, not necessarily 
true that this value will exceed the one for normal two-dimensional flow. 

This rather crude description the secondary flows may help explain 
some the essential features theflow mechanics, but more 
adequate description can vortex generation. Some con- 
cepts bearing this approach are described Hawthorne (6) and Squire 
and Winter (10) but more study needed developthem for the 
problem this investigation. 


EXPERIMENTAL INVESTIGATION 


Experiments were variation the local drag coefficient 
along cylinder. diam hard copper tube was mounted vertically the 
center 12-in. wide glass-walled flume which flow water was main- 
0.8 ft. Theflow adjusted togive Reynolds Numbers 
mid-depthin the range between 18,500 and 22,500,a region normal 
two-dimensional drag constant. With the Reynolds Number inthis 
range, and the form and held constant, these factors were eliminated 
source variation the measured drag coefficients. Fig. schematic 
diagram the test arrangement. 

Since the flow had free surface might appear that the drag coefficients 
would influenced the free surface and hence the Froude Number should 
considered. This would course true the total drag coefficient were 
under consideration because the drag produced the wave the surface could 
very important part the total drag. The investigation reported here 
deals with the local drag coefficient determined from pressure measurements 
the part the cylinder that was completely submerged. believed that 
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the free surface had negligible effect onthe pressures measuredat the cylin- 
drical surface except for points near thefree surface. is, therefore, believed 
that the free surface had effect the results reported here. 

Apparatus Control the velocity distribution was variable 
shutter gate. Variation the gate opening from top bottom permitted the 
resistance varied from maximum near the bottom minimum the 
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FIG, 1.—DIAGRAM TEST ARRANGEMENT 


surface. The effectiveness the gate producing linear velocity variation 
was enhanced the inclusion serrated screen. 

The cylinder was mounted aluminum plate suspended four thin alu- 
minum frame. (Fig. The lower plate was free deflect 
force was applied the cylinder. Deflection the lower plate introduced 
flexural stresses the aluminum straps, the strain from which was picked 
system electrical strain gages, thereby enabling total force measure- 
ments obtained. 
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FIG, 2.—CYLINDER MOUNTING ASSEMBLY 
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The cylinder was attached the lower plate sleeve fitted into collar 
fixed the lower plate. was possible rotate the cylinder through 360° 
about its longitudinal axis. Located along vertical line the cylinder were 
piezometers, each which was connected internally toa 1/8 in. copper 
tube. The tubes were brought out through the top the cylinder andconnected 
manometer board from which the local pressure was measured for each 
piezometer. Fig. shows the cylinder and the amounting assembly. 

Experimental procedure.—The experimental procedure involved measure- 
ments the velocity distribution, local pressure distribution, and total drag 
force. The longitudinal velocity profiles were measured with Prandtl pitot 
tube. Several transverse profiles were also obtained insure that side wall 
effects were negligible. Pressure distributions around the cylinder were ob- 
tained with the piezometers rotating the cylinder through 360° increments 
Total drag force measurements were taken several times during each 
run. All these measurements were taken four different values the ve- 
locity gradient parameter. Total drag measurements were made check 
the integrated pressure measurements. The measured total drag included, 
course, the wave drag the free surface. approximate evaluation the 
wave drag was made, based the profile the wave surface around the cyl- 
inder. The sum the estimated wave drag and the pressure drag obtained 
from integration the pressure measurements agreed satisfactorily with the 
total drag force measured the electrical strain gages. 

The evaluation the local drag coefficient various points along the cylinder 
was determined from the pressure distributions and the measured velocity pro- 
files. the method determining the local drag coefficient, consider 
section the cylinder. The pressure distribution determined from the 
piezometer the integrated tofind theforce the particular sec- 
tion. the functional relation betweenthe pressure and the angle couldnot 
expressed Eq. 4was evaluated graphical integration. Typical 
the curves used for this integration Fig. This plot includes 
curve for each piezometer. 

Having completed the integration the pressure distribution section, 
the corresponding local velocity was taken from the plot the velocity profile, 
Fig. The local drag coefficients were computed various positions along 


the cylinder for different values the velocity gradient parameter inthe man- 
ner described above. 


RESULTS INVESTIGATION 


Local drag coefficient.— The velocity gradient very definitely affects the local 
drag coefficient. The exists may best seen Fig. Here the 
local drag coefficient plotted against the distance from bottom feet for 
different values the gradient parameter. the velocity distributions were 
uniform, there would the local drag coefficient the velocity 
and pressure distribution would constant along the cylinder. This curve in- 
dicates that the local drag coefficient does vary with position and also with the 
value the slope the velocity profile. 

Consider first the variation with position along the cylinder. Fig. indicates 
that drag maximum the bottom the cylinder and 
minimum With velocity gradient, this corresponds 
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LOCAL DRAG COEFFICIENT 


PER. SEC. 


FIG, 


FIG, 4.—VELOCITY PROFILES 
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high coefficient the low-velocity end and low coefficient the high- 
velocity end. This variation coefficient was,as pre- 
viously described, terms the secondary flows. 

There is, also, variation the local drag coefficient with the steepness 
the velocity gradient. From the curves Fig. can seen that the 
slope parameter increases from 0.058 the variation the local drag co- 
efficient becomes more pronounced. This the fact that the greater 
the velocity gradient, the greater the pressure variation the wake zone, 
and the greater the upward flow the wake zone. 

illustrate the upward flow that exists the wake, dye was injected 
through one the copper tubes issuefrom piezometer located 
the wake zone. The result this qualitative demonstration presented 
Fig. The dye jet may seen shooting out from the piezometer for dis- 
tance about 1/4 the cylinder diameter and then turning upward where 
hugs the back the cylinder. Although there some subsequent diffusion 


* 


FIG, 6.—UPWARD FLOW WAKE, 


the dye into the surrounding fluid swept downstream, the existence 
the upward flow emphasized the fact that only trace the dye appears 
below the piezometer opening for the downstream region. 

Average drag coefficient.—In general, may concluded that the average 
drag coefficient for cylinder witha longitudinal velocity gradient less than 
the conventional value 1.2 for the chosen range Reynolds Numbers. The 
primary reason for this that not two-dimensional. The presence 
the longitudinal flow the wake zone very definitely indicates three-dimensional 
flow characteristics. 


The values the average drag coefficient agree quite well with the approx- 
imate values given for cylinders finite length. For values 20, the 


average drag coefficient varies from 0.74 0.90. (3) The results this in- 
vestigation more nearly approach these values the steep velocity gradients, 
however, the causes are not directly related. 


SIGNIFICANCE FOR ENGINEERING PROBLEMS 


The investigation has shown that the longitudinal velocity gradient can effect 
the value the drag coefficient and that longitudinal flow induced the 


~ 
j 
= 
ihe, 


wake the cylinder which appears related the effect. 
compare the values the velocity gradient parameter used the investigation 
with values that might expected some engineering problems. the ex- 
perimental investigation the parameter AV/AY D/V had the following values: 
0.117, 0.077, 0.061, and 0.058. 

For comparison, some values the velocity gradient parameter were cal- 
culated for waves based laboratory measurements particle velocities 
waves Marlowe. (4) assumed thata pier diameter 0.05 times 
the water depth, the wave length 2.5 times the water depth, the double ampli- 
tude 0.11 times the water depth, and the velocity gradient mid-depth 
used, the gradient parameter will 0.10. Use the velocity gradient 0.8 
depth and all other factors the same gives gradient parameter 0.135. Ina 
wave greater height, that with double amplitude 0.145 times the water depth, 
and using the velocity gradient 0.8 depth the gradient parameter 0.25. 

If, for simplicity, one assumes that the wind velocity distribution near the 
ground may the 1/6 power the elevation, (1) approximate 
value the gradient parameter may computed fora stack. one considers 
diam stack 100 high, and bases the gradient parameter the gradient 
mid-height the gradient parameter will 0.05. one considers the lower 
portion where the gradient steeper, for example, the lower ft, the gradient 
parameter would 0.125. These values indicate that the gradient parameters 


used the investigation are inthe significant for some engineering 
problems. 


CONCLUSIONS 


this exploratory investigation the velocity gradient effect, test condi- 
tions were simplified the following ways: 


The ratio cylinder diameter water depth was essentially constant. 

The Reynolds number range was chosen that the conventional drag 
coefficient was essentially constant. 

The velocity variation was approximately linear giving constant value 
for the gradient all points along the cylinder. 


Based the results obtained under these restrictions, the following con- 
clusions appear justified: 


longitudinal velocity gradient the flow approaching circular cylin- 
der induces longitudinal currents the wake zone and along the upstream ele- 
ment, thus, altering the local drag coefficient various positions along the 
cylinder. 

The local drag coefficient reduced sections subject highest 
stream velocity due the rise pressure fluid supplied into the wake 
from the regions lower stream velocity. 


The magnitude the above effect increases the magnitude the ve- 
locity gradient increases. 


These exploratory studies have indicated that the velocity gradient effect 
significant number problems engineering importance, and has 
opened several questions requiring further investigation. Further study 
anticipated evaluate the influence some the factors that were excluded 
order simplify this exploratory investigation. 
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UNSTEADY FLOW GROUND WATER INTO SURFACE 


William Haushild! and Gordon Kruse” 


SYNOPSIS 


Prediction the water table position and the amount water discharged 
where ground water flowing from aquifer surface reservoir has not 
been exact. The nonlinear partial differential equation that describes the shape 
the water table difficult solve. 

Approximate solutions obtained two different methods for the nonlinear 
equation are presented. Both approximate solutions agree better with experi- 


mental results than does the exact solution the simplified linear equation 
for the flow. 


INTRODUCTION 


The solution one-dimensional boundary problems ground water flow 
has many applications. For example, the increment water added bank 
storage when surface reservoir filled, taken from bank storage when 
reservoir emptied, can considerable and the amount important ina 
determination the reservoir hydrology. The same solution also applies 


Note.—Discussion open until December 1960. extend the closing date one 
month, written request must filed with the Executive Secretary, ASCE. This paper 
part the copyrighted Journal the Hydraulics Division, Proceedings the Ameri- 
can Society Civil Engineers, Vol. 86, No. July, 1960. 

Joint contribution from Quality Water Branch, Geol. Survey, Dept. 
Interior and Soil and Water Conservation Research Div., Agric. Research Service, 
Dept. Agric. 

Hyd. Engr., Quality Water Branch, Geol. Survey, Stationed Colorado State 
Univ., Fort Collins, Colo. 

Agric. Engr., Soil and Water Conservation Research Div., Agric. Research Service, 
Dept. Agric., Stationed Colorado State Univ., Fort Collins, Colo. 


| 
Vi. 


certain unsteady cases ground water flow caused the digging drains 
and bank storage adjacent streams that undergo changes water level. 
This paper presents four approximate solutions for the problem. Solutions 
and are first approximations which the effect the unwatering the 
aquifer the area available for the flow ground water neglected. Solu- 
tions and are second approximations that take the effect unwatering 
into account. 

The solutions presented inthis paper are based upon the following assump- 
tions: 

The saturated aquifer infinite extent and overlays impermeable 
layer zero slope. 

The Dupuit-Forchheimer assumptions hold. 

The change reservoir water level instantaneous. 

The aquifer composed isotropic and homogeneous material. 


SOLUTION 


Consider the flow into reservoir shown Fig. The flow through 
cross section width the distance from the origin 


which the permeability the aquifer. 
let 
(2) 
then 
The equation continuity for the strip width and time increment 
which the specific yield the aquifer. 
Or, use Eq. 
32h 


can discarded from Eq. The continuity equation (Eq. then gives the 
following linear differential equation: 
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Eq. solved Carslaw and Jaeger3 applies the flow heat 
solids. applied Eq. the case constantly discharging surface 
drains. solved equation similar Eq. numerical methods. 


Land Surface 


Water Table 


Reservoir 


FIG, 1.—GROUND WATER DISCHARGE INTO BODY SURFACE WATER 


The results obtained concurred with data obtained from Hele-Shaw viscous 
flow model. 
The solution Eq. given Carslaw and Jaeger 


which the error function, 


The integral tabulated 

This solution When the drawdown small compared the thickness 
the aquifer, solution approximate solution Eq. satisfies the 
boundary and initial conditions: 


“Conduction Heat Solids,” Carslaw and Jaeger, First Edition, 
Oxford Univ. Press, 1946, pp. 50-51. 

Quantitative Method for Determining Ground-Water Characteristics for Drain- 
age Design,” John Ferris, Agricultural Engineering, 31:6, 1950, 285. 

“Unsteady Flow Porous Media Means Hele-Shaw Viscous Fluid Model,” 
David Todd, Transactions, Amer. Geophysical Union, 35:6, December, 1954, 

“Tables the Error Function and its Derivative,” Natl. Bur. Standards Applied 
Math. Series No. 41, October, 1954, 
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SOLUTION 


second approximation the true solution Eq. may obtained 
applying the method This method uses each successive approxi- 
mation approach the correct solution. obtain solution the nonlinear 
terms Eq. are computed from the first approximate solution, Eq. and 
substituted into Eq. known functions. The differential equation, thus 
modified, then solved, again being subject the initial and boundary condi- 
tions. The process solution requires that particular integrals found for 
the known functions. 


particular integral for the term 
Similarly particular integral po, for the term 
The expanded form this particular integral 
4at 


approximate solution the modified form differential Eq. that satisfies 
the initial and boundary conditions 


This solution II, inwhich the first approximation, Eq. altered cor- 
rect for the initial and boundary conditions. 


plot Eq. shown Fig. with the dimensionless parameter h/H 
that the solution compared with experimental data. The solution for the 
linear case, Eq. also shown Fig. 

Keller and obtained measurements the drawdown curve ina 
laboratory experiment sand filled Water was drained from the sand 


“Memoire sur theorie des equations aux derivees partielles method des 
approximations successives,” Emile Picard, Liouville’s Jour. Methematiques, 
Series Vol. 1890, pp. 145-210. 

“Model Study Interceptor Drains, Jack Keller and Robinson, Proceed- 
ings, ASCE, September, 1959, pp. 25-40. 
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All solutions are computed for 
reservoir level lowered 
tom the aquifer. 


FIG, 2.—APPROXIMATE SOLUTIONS GROUND WATER 
EQUATIONS AND EXPERIMENTAL DATA 


ia 
| 
A 
/ 
ie 
ae 
a¥- 
-04 
| 
4 
f 


July, 1960 


through 4-in. perforated pipe located the floor the flume. The entire 
flow water was intercepted the drain. These conditions simulated flow 
into surface reservoir where the water drawn down the bottom the 
aquifer 0). The shape the water table was determined before the draw- 
down curve reached the upstream boundary the flume. Results were, there- 
fore, the same for aquifer infinite extent. The sand used had values 
0.034 fps, and 25.5% and was practically homogeneous. The data 
Keller and Robinson are compared with the theoretical approximation, Eq. 13, 
Fig. 


SOLUTION 


alternative second approximation for the bank storage case may ob- 
tained method equating flows. first approximation obtained 
solution and the unit discharge, computed from it. This flow then 
substituted into Eq. and improved solution obtained integrating the 
resulting expression. The first approximation will solution III. ob- 
tained solving the linearized differential equation 


subject the 


The significance the notation used this case shown Fig. 
The required first approximation, solution III, then 


Solutions and III satisfy the same equations and boundary conditions with only 
the reference point changed. These solutions, therefore, are identical. 


Water Table 


Reservoir 


FIG, 3.—GROUND WATER DISCHARGE INTO BODY SURFACE WATER 
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SOLUTION 


first approximation the flow through cross section unit width 


which D=d and given Eq. 15. 
the drawdown taken into account, the expression for flow 


The constant may evaluated for the boundary condition that when 


Then, substitution 


Solution obtained from this relation the form 


from the same datum that they may compared directly, Eq. be- 


comes 


For the case complete drawdown takes the form 


This equation shown Fig. for comparison with experimental data. 

comparison these solutions shows that the second approximations 
and show closer agreement with the experimental data than the first ap- 
proximations and Fortunately, the best agreement obtained from so- 
lution IV, which the easiest use. The comparison made here with the 
extreme case complete drawdown. 

All the solutions are approximations based the assumption that very 
large compared This assumption least valid for the case complete 
drawdown. Therefore, the results from using this solution with partial draw- 
downs should good. 


| 


Solutions through were derived for the case inflow into reservoir 
upon lowering the reservoir water level. However, the same equations and 
solutions the case outflow considered negative for flow 
into the aquifer. 


COMPUTATION 
Curves similar those Fig. can used compute for given 


aquifer values andt have been measured the field. From 
these values H/D can computed and curves drawn for solution solution 


TABLE 1.—ERROR DETERMINATION 


h/H Solution Solution IV, 
percent percent 

-0.12 

+0.06 


using h/H ordinate and abcissa. For given h/H, can 
computed from the measured value and the value obtained 
from the curve. 

Solutions and were used todetermine bythe above method for several 
the points from Keller and Robinson’s data.8 These values were compared 
with the value computed from laboratory measurements and 
Table indicates the percentage error determined these two solu- 
tions. 
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DISCUSSION 
DISCHARGE FORMULA FOR STRAIGHT ALLUVIAL 


Discussion Garde and John Bogardi 


GARDE, A.M. authors are complimented for their 
work concerning the discharge formula for the straight alluvial channels. They 
have shown that the mean velocity flow straight alluvial channel 
given the equation 


then C,, and are not constants but they depend the size the bed ma- 
terial and the regime flow. According the analysis carried Messrs. 
Liu and Hwang, can vary between and 287, can vary between 0.35 and 
0.71, and can vary between 0.30 and 0.57. For the large volume flume 
data, the correlation obtained seems remarkable. 

During 1956-59, the writer carried out certain concerning 
the velocity distribution and the sediment transport alluvial channels. One 
the significant findings these investigations was the fact that velocity dis- 
tribution and sediment transport are very much dependent whether the bed 
the channel plane, has ripples and dunes The work Messrs. 
Liu and Hwang confirms this finding. 

order find out whether this new formula proposed Messrs. Liu and 
Hwang can applied natural channels and rivers, the writer has carried 
out certain preliminary analyses and the results are reported herein. The data 
used for this analysis cover wide ranges the discharge the mean velocity 
the hydraulic radius and the bed material size Table gives sum- 
mary these data. 

The procedure followed was: 

determine the regime flow according Criterion. 


Note.—This paper part the copyrighted Journal the Hydraulics Division, 
Proceedings the American Society Civil Engineers, Vol. 86, No. July, 1960. 

Civ. Engrg., Univ. Roorkee, Roorkee (U.P.), 

“Total Sediment Transport Alluvial Channels,” Garde, Ph. Disserta- 
tion submitted Colorado State Univ., Fort Collins, Colo., January, 1959. 
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compute the mean velocity flow the use Eq. 


was found that several points from these field data fell the region 
‘Transition’ Fig. (given the authors). When this predicted regime was 
used and the corresponding values Cg, and were chosen, the values 
the mean velocity could computed. was found: 


For points falling the dune regime, the observed velocity, most the 
time, was greater than the computed velocity. 


For the points falling intransition regime the observed velocity was less 
than the calculated velocity. 


Considering, for example, the Punjab Canal’s data, points out fell 
the transition regime and for these points the computed velocities were 


TABLE DATA 


(1) (3) (4) (5) (6) 
3.3 
Sind Canals 4-12 300- 
3.4 9000 
Ganga Canal (India) 3.4- 7400- 
12.4 4.0 9700 
Indus River 17.8- 6.8- 3,35000- 
26.8 7.0 5,80700 
Tiber River 2.0- 4500- 
3.1 11600 
Rakaia (New 5.0 4500 
Opihi (New 1.1 
Wanganui (New 1.4 1630 
2.3 
Canals 2.3 
2.6 


and Design Stable Channels Alluvial Materials,” Simons, 
Ph. Dissertation submitted Colorado State Univ., Fort Collins, Colo., May, 1957. 
“Slope Discharge Formulae for Alluvial Streams and Rivers,” E.C. Schnacken- 
berg, Proceedings, New Zealand Inst. Civ. Engrs., Vol. 37, 1951. 
Agric. SCS-TP 55, August, 


greater than the observed velocities 35% 125%. For the points that 
fell the dune regime, the computed velocities were always smaller than the 


observed velocities 14% 43%. Similar tendency was noticed the case 
Sind Canal’s data. 
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DISCUSSION 


One the possible reasons for such (in the case data fal- 


ling ‘transition’ regime) could that the regime predicted 

criterion not correct. fact, the authors have stated that this criterion 

based the flume data for which the Froude numbers are much larger than 

those for natural streams. The accuracy the results obtained the Liu- 

Hwang formula depends the correct prediction the regime and therefore 

the criterion for determining regime flow very important factor. 

pointed out the authors, criteria predict the 
regimes flow with fair degree accuracy for the flume data. However, 
when applied the natural river data with larger depths and flatter slopes, 
the criteria encounter serious the case natural river data, 
because high values shear, attains large value and thus the points 
fall transition antidune region. Actually these points belong ripple 
dune region. This discrepency occurs because the Froude number not taken 
into consideration while determining the regime. 

From the point view dimensional considerations can shown that, 


asa approximation, the ratio heightof dunes length dunes depends 


that the Froude number important factor the determination regimes 
flow alluvial channels. Similarly, important the determi- 
nation velocity distribution alluvial For these reasons one 
would expect that and will significant parameters for establish- 


ing the criterion for the determination regimes flow. This was done 
Garde and and this criterion seems quite satis- 
factory. applicable flume data well natural river data. 

Since the Froude number, which one the parameters determining re- 


gimes flow, contains the mean velocity flow, this criterion 


cannot used while determining the mean velocity flow. However, this 
criterion can used with advantage the present case where the purpose 
study the applicability equation natural rivers and canals. 

For these the regimes flow all the data were determined ac- 


cording criterion; several points which fell the ‘transition’ 


regime according criterion, now fell the ripple-dune regime. 


This made considerable difference the computed mean velocity flow 
according Table shows the comparison between the observed mean 


velocity and the computed mean velocity for the regime predicted the two 
criteria. 


“Hydraulic Criteria for Sand Waves,” Langbein, Transactions, Amer. Geo- 
physical Union, 1942, pp. 615-618. 

Matunashi, Personal letter. 

“Characteristics Bed Forms and Regimes Flow Alluvial Channels,” 


Garde and Albertson, Civ. Engrg. Sect., Report CER RJG Colorado State 
Univ., Fort Collins, Colo. 
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TABLE B.—COMPARISON BETWEEN OBSERVED AND COMPUTED MEAN VELOCITY 


Canalorriver Calculated 
Transition regime 
(1) 


Sind Canals 


Ripple-dune regime 


(4) 


Tiber River 


| PUNJAB CANALS 
INDUS RIVER 
TIBER RIVER 

RAKAIA RIVER 
OPIHI RIVER 
WANGANU! RIVER 
GANGA CANAL 
MANGLES RIVER 
SIND CANALS 
MOUNTAIN CREEK 
U.S. CANALS 


LINE PERFECT 
AGREEMENT 


COMPUTED VELOcITy FPS 


OBSERVED VELOCITy FPS 


FIG, OBSERVED AND COMPUTED MEAN VELOCITIES, 
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Thus can seen that computed velocities according ripple-dune are 
better agreement with observed velocities. Fig. shows the comparison 
between observed mean velocity flow and computed mean velocity for the 
various data. can seen that the general tendency for the observed ve- 
locity greater than the computed mean velocity. From Fig. can aiso 
seen that 35% the data have less than 25% error. 

slight adjustment the curve better agreement could ob- 
tained. The following procedure was adopted. was assumed that the values 
and obtained from the authors’ graphs were correct and that adjustment 
could made for values such manner that mean velocity 
flow equal the calculated mean velocity flow. Thus 


New values were computed this manner for all the data available 
the writer. When there was more than one observation for the same size 
bed material and the same regime flow, the average new value was 
found. Fig. shows the variation with for the ripple-dune regime. 
Lane and Carson’s data used Messrs. Liu and Hwang are also plotted 
Fig. Tentatively mean curve passing through these points drawn. 
the size bed material 0.5 0.6 mm, this new curve very close 
the one proposed the authors. However, for coarser material the proposed 
curve gives lower values than the one given the authors’ curve. More 
data are needed plotted this figure, especially for coarser material, 
order confirm the location the proposed curve. 


JOHN attempt the authors investigate this ex- 
tremely involved problem very much appreciated. The writer agrees with 
the statement that with the effect sediment transport upon resistance being 
unknown, empirical relationships rather than theoretical solutions should pro- 
vide the correct approach. The soundness this attitude corroborated 
the successful study the authors. 

almost complete relationship between the three invariants and 


these two figures, the points representing identical values obviously 


always result from the same experiment, which the particle size, the spe- 
cific gravity the sediment, and the water temperature are uniform, 
which there are but slight changes the water temperature only. Each the 
parallel straight lines Figs. and sloping 0.555 and 0.565, 
respectively, represent the result experiment presumably carried out 
under different conditions observation. Nevertheless, positive relationship 
can established Figs. and between the coefficient “A” and the invari- 


ant which corroborates the general validity the relationship, termed 


empirical correlation the authors. 


Head, Hydr. Lab. Research Inst. for Water Resources Development, Budapest, 
Hungary. 
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DISCUSSION 


the writer’s opinion the explanation given the authors asto the compo- 
sition the invariant group slightly ambiguous. According this, the 
composition should supported the following considerations: 


conform with the existing knowledge boundary resistance. 
correlate the data consistently. 
Namely all the variables 


V. =,/g R,S, 


5,—-, F= 


with the exception the unknown velocity are contained the invariants 

principle, there could objection using variable, respectively 
invariant, differing essentially from Subsequently the authors themselves 


admit that the good applicability the invariant group (namely that the 


points representing identical values should straight lines) has been at- 


tained the appropriate selection values and 
fact, wide variety solutions can conceived for function the 
form 


The circumstance, that the above invariants will involved the function, 
has been illustrated the authors themselves the aid 
ysis. the basis few measured data the writer has thus temporarily de- 
termined the function Eq. with fair accuracy six-variable correlation. 
appears that with more data available the value respectively 
could determined graphically with very great accuracy. 

Regardless these considerations, the writer wishes state that the au- 
thors have made very significant advances this difficult problem. The writer 
agrees with the authors’ statement concerning the exponents and his 
opinion the method used for determining and such already contain 
the effect various regimes movement. likely, therefore, that and 
are related the invariants that are characteristic for movement. 

The attitude the authors towards the application dimensionless analysis 
may give rise misunderstanding. Actually, the use dimensional analysis 
alone does not guarantee that there has not been any important variable omitted 
from the investigation. 
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connection with the authors’ comments its original and pres- 
ent form, the writer would like point out, that representing the relation- 
ship 


the tractive force can all probability determined even simpler form. 
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THE SETTLING PROPERTIES 


Discussion Charles Gunnerson 


CHARLES GUNNERSON, ASCE.--The author has pointed out that the 
analyses settling velocities various sewage effluents shown his Fig. 
were made connection with studies marine disposal the vicinity Los 
Angeles. During the same period observations were made the behavior 
these effluents the receiving waters. The results the marine studies have 
been previously reported.2 sufficient here summarize the data relating 
the sedimentation effluent suspended solids sea water. 

The sanitary significance sedimentation suspended solids from the ef- 
fluent the Hyperion Treatment Plant lies the concurrent reduction en- 
teric bacteria populations surface waters. Accordingly, observations were 
made the rate disappearance coliform bacteria from the surface. The 
probable effects dilution and mortality were determined and the effect 
sedimentation was computed. 

was found that, the average, coliform bacteria were removed from 


surface sedimentation rates such that 90% removals occurred 
follows: 


Hyperion secondary effluent (plus elutriation overflow)—21.0 hr. 
Hyperion primary effluent (plus elutriation overflow)—5.3 hr. 
Orange County Sanitation Districts primary effluent—2.0 2.4 hr. range 


given because some uncertainty indetermining the effect di- 
lution). 


These results are based upon samples taken from the upper 1-ft layer the 
water over periods day. 

order compare the field results with the author’s settling tube studies, 
the data his Fig. have been replotted show suspended solids reductions 
assumed that the sedimentation velocities represent removals from slice 
the depth. further assumed that such removals can compared, 


least gross basis, with removals from the essentially homogeneous 
upper 1-ft layer the ocean. 


December, 1959, Ronald McLaughlin Jr. 

Civ. Engr., Bur. Sanitation, City Los Angeles, 

“Sewage Disposal Santa Monica Bay, California,” Gunnerson, Proceed- 
ings, ASCE, Vol. 84, No. February, 1958, 1-28, and Transactions, ASCE, Vol. 
124, 1959, pp. 823-851. 

“Studies Coliform Bacteria Discharged from theHyperion Outfall, Final Bacter- 
iological Report,” Rittenberg, Allan Hancock Foundation, Univ. Southern 
California, Los Angeles, Calif., 1956. 
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Fig. shows both the laboratory and field data for the removal suspended 
solids with time for Hyperion and Orange County Sanitation Districts effluents. 
Although there are minor differences the groups effluents tested, the 
rates sedimentation the ocean are clearly greater inthe settling 
tubes. Qualitatively, least, the effect flocculation the colloidal fraction 


the suspended solids appears account for the greater part solids re- 
moval the ocean after the first few hours. 


uJ 
| | 


TIME, SECONDS 


FIG, SEDIMENTATION RATES 
SETTLING TUBES AND THE OCEAN, 


The flocculation stream-carried suspended solids discharged into salt 
water has been studied marine geologists for many years. 
considers that the variable composition source rocks and diagenetic changes 
after deposition the ocean floor are probably more important factors inthe 
distribution clay minerals marine sediments off the California coast than 


«The Sea Off Southern California,” Emery, Wiley, 1960. 
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DISCUSSION 


are differential flocculation and deposition. Nevertheless, van Andel and 
explaining the near- and off-shore distribution illite and montmorillonite, 
respectively, the Gulf Paria. 

Some the aspects Whitehouse’s results studies differential floc- 
culation and settling clay minerals appear applicable sedimentation 
effluent solids the ocean. found that the median settling velocities in- 
creased about two orders magnitude upon flocculation saline waters 
with parts per thousand chlorinity for the various clays tested. Ini- 
tial settling rates increase and 22% for kaolinite and illite, respec- 
tively, the chlorinity increased from 0.5 about and are unaf- 
fected further chlorinity. contrast, the chlorinity in- 
creases from 0.5 initial settling rates for montmorillonite are in- 
However, the velocities are still the order 1/10 those the other clays. 
The effect temperature reduction from 21° (this about the sea- 
sonal range Santa Monica Bay the upper 200 water) decrease 
settling rates the three clays about 30%. Ions heavy metals increase 
the settling rates; the highest increase 85% was reported for montmorillo- 
nite the presence 10-4 molal The presence carbohydrates 
increases the settling velocity montmorillonite 35% and 25% 0.5 and 
chlorinity, respectively. Humic acids proteins cause decreases 
clay concentration decreases, there small but significant increase set- 
tling velocity. Other factors which influence settling are pH, clay concentra- 
tion gradients, and the initial state division the clay. 

The fate dispersed organic solids not well characterized. Some 
the organics are undoubtedly absorbed onto the inorganic particles, some may 
flocculated and settled, while the remainder may remain dispersed for 
considerable time. Fox, have found colloidal otherwise 
dispersed organic detritus mid- and deep-ocean waters well near- 
shore waters. follows that the persistence this materialin the water col- 
umn determined the rate bacterial decomposition and mineralization 
rather than flocculation and settling. 


“Recent Sediments the Gulf Paria,” Tj. van Andel and Postma, Reports 
the Orinoco Shelf Expedition, Vol. North Holland Pub. Co., Amsterdam, 1954, 

Sediments and Their Environments, AP1 Research Project 51, Reference 52-51, Scripps 
Inst. Oceanography, 1952, pp. 23-32. 

53-43, Scripps Inst. Oceanography, 1953, pp. 31-38. 

“Diagenetic Modification Clay Mineral Types Artificial Sea Water,” 
Whitehouse and McCarter, Clays and Clay Minerals, Natl. Acad. Sciences- 
Natl. Research Council, Washington, C., Publication 566-1958, pp. 81-119 (Includes 
reference and above). 

“Differential Settling Tendencies Clay Minerals Saline Waters,” 
Whitehouse, Jeffrey, and Debbrecht, Proceedings, Seventh Natl. Conf. 
Clays and Clay Minerals, Pergamon Press, 1960, pp. 1-79. 

“Marine Leptopel, its Recovery, Measurement, and Distribution,” Fox, 
Isaacs, and Corcoran, Journal Marine Research, No. 1952, pp. 


“particulate Organic Detritus,” Fox, Geol. Soc. America. Memoir 67, 
Vol. 1957, pp. 383-390, 
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The previously cited works are not directly applicable waste disposal. 
Nevertheless, reasonable expect that least portion the effluent 
solids will respond various environmental factors manners similar 
those natural compounds. Certainly, differential flocculation phenomena 
must considered evaluating the rate and location deposition effluent 
suspended solids and their constituents. The converse this proposition, that 
dispersed, unflocculated solids transported for greater distances may 
have similar significance. 

One the conclusions derived from the marine studies the disappearance 
rates coliform bacteria from Hyperion discharges was that primary effluent 
was more amenable final treatment the ocean than was secondary ef- 
fluent.2 Two approaches have been used attempting verify these results 
under controlled laboratory conditions. 

the first, various sea water dilutions the effluents were placed large 
cylindrical columns, and tests were made over period time coliform 
densitites different depths. was hoped that the tests would indicate rate 
sedimentation the bacteria and consequent reduction bacterial popu- 
lations the surface. such reduction was found, probably because con- 
vection currents being set the column bysmall temperature gradients 
the room and because the motility the bacteria themselves. 

The second series tests involved centrifuging various sea water dilutions. 
(The writer indebted Mysels the University Southern California 
for suggesting this technique). The tests were made room temperature with 
the centrifuge operating 1000 g’s for min. Qualitatively, primary effluent 
mixtures showed greater reductions coliforms than did secondary, thus con- 
firming the receiving water studies. was also found that greater reductions 
were obtained with higher dilutions; this conforms with the work Whitehouse, 
al., mentioned earlier. 

was possible determine the suspended solids reduction with assurance 
only for the 20-to-1 dilution primary effluent since the absolute solids con- 
centrations the other dilutions were too small accurately measured 
with standard methods. The solids reduction was about two-thirds that the 
bacterial reduction. 

now have three converging sets data, the marine studies, the cen- 
trifuging experiments, and the settling tube investigations reported Mr. 

McLaughlin. The differential settling characteristics various sewage ef- 
fluents are clearly demonstrated, particularly when effects flocculation 
the effluent solids sea water are included. 

Associated with the utility centrifuging sea water, effluent mixtures 
are also indicated. However, obvious that this technique requires good 
deal experimental and analytical study before quantitative results can ex- 
pected that will have the same the very expensive receiving water 
studies. any event, appears desirable make laboratory tests settling 
characteristics effluent solids which may form bases design for sewage 
treatment plants. Such experiments may well designed include the floc- 
culating other effects the actual receiving water, they fresh water 
salt water. 
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EARLY HISTORY HYDROMETRY THE UNITED 


Discussion Murray Blanchard, Arthur Frazier, and Stevens 


MURRAY BLANCHARD, ASCE.—The author has compiled avery inter- 
esting history hydrometry. Under the heading Great Lakes Investigations 
has given complete resume the equipment and methods developed and 
used the engineers the United States Lake Survey this work. 

The writer had opportunity make use those pioneer developments 
the following hydraulic investigations the Great Lakes: 


Discharge from Lake Superior the St. Marys Sault Ste. Marie 
Mich., for the United States Lake Survey. 

Dischargefrom Lake Huron the St. Clair River Port Huron, Mich., 
for the United States Lake Survey. 

Discharge the Detroit River between Lakes St. Clair and Erie 
Detroit, Mich., for the United States Lake Survey. 

Discharge the St. Lawrence River the foot Lake St. Francis near 
Valley Field Ontario, Canada, for the White Company New York 
connection with Water Power development investigations. 

Discharge the Chicago Sanitary District’s Main Channel from Lake 
Michigan Lemont, for The Sanitary District Chicago connection 
with United States Lake Level 


the hydrometric investigations the United States Lake Survey referred 
the author there should added discharge measurements the outflow 
Lake Superior made and 1902 Johnson and Edward Wilson 
from the International Bridge Sault Ste. Marie, Mich. 


ARTHUR finds hard realize that this paper 
American engineering history has been written person who has lived 
America for only twelve years. doubt most the native-born Americans 
who read found much information therein that was entirely new them. 
accordingly with humble spirit that this discussion offered. 

only one detail find reason disagree with data contained the pa- 
per. The author states that the first discharge measurement made State 
Maine was made Hamlin Maine State College who measured the 
Penobscot River 1884, whereas George Vose the Massachusetts Insti- 
tute Technology 1885 that: 


January, 1960, Steponas 

Retired Civ. Engr., Libertyville, 

Blanchard, Hydraulics Chicago Sanitary District’s Main Channel, Journal 
the Western Soc. Engrs., Vol. 25, No. 13, Sept. 1920, pp. 471-524. Discussion: 
Vol. 26, No. August, 1921, pp. 300-303, Chicago. 

Retired member the Geodetic Survey. 


Sketch the Life and Works Loammi Baldwin,” George Vose, Boston, 
1885. 
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“In 1835, several prominent gentlemen Maine became interested 
the development the water power the Androscoggin River 
Brunswick; and Mr. Baldwinwas employed make the necessary meas- 
urements and computations .... From the several gaugings, 4,000 cu- 
bic feet per second were reckoned safe the discharge which might 
relied upon constant through the year for mill power; and with forty 
feet fall the river this place, Mr. Baldwin concluded that the Andro- 


scoggin Brunswick offered power unsurpassed any New England 
that time occupied.” 


matter fact, seems quite unlikely that even Baldwin’s measurement 
was the first one made Maine. contemporary, and doubt close friend 
his, Ithamar Beard, made surveys for the erection large cotton fac- 
tory Saco, Maine, 1831. Not long afterward Beard advertised (in the 
Lowell Courier during October, 1840) that he, civil engineer, offered his 
services “in surveying measuring and computing the quantity 
water supplied streams ..., and generally, any business pertaining 
engineering.” His records streamflow measurements, and those similar 
engineers his time, would provide interesting comparisons with present-day 
data they could found. more than likely that many Beard’s meas- 
urements would have been made earlier than any presently record for the 
states which lived (Maine, New Hampshire, and Massachusetts). 

writing history such this, one always faced with the problem 
what items might omitted for the sake brevity. That circumstance al- 
ways provides almost unlimited area for discussion purposes. Advantage 
will taken thereof the comments that follow. 

Perhaps Nicolas Nicollet was singled out among the explorers the 
Upper Mississippi River because the magnificent map made the area 
that early data was not, however, the first government- 
sponsored explorer have undertaken the job. Colonel Zebulan Pike was or- 
dered search for the head-waters that river 1805. General Cass and 
later (1832) Henry Rowe Schoolcraft also were sent that mission. those 
three predecessors Nicollet, only one, Henry Rowe Schoolcraft, actually 
reached Lake Itaska (the nominal “headwaters” the Mississippi). re- 
markable coincidence that still earlier explorer this area also bore the 
name Jean Nocolet. 1634 this earlier Jean Nicolet discovered Lake 
Michigan, and even undertook journey the Fox River Wisconsin 
quest the “Father Waters,” but did not far, probably only far 
Lake Winnebago, and then turned back, only day’s journey away from his goal. 
For the sake the record, however, question might raised regarding the 
given name attributed “Jean” Nicholas Nicollet the author. Most Ameri- 
can sources have translated his first name “Joseph,” rather than “Jean.” 

From historical viewpoint would seem that report made April 
1808 canals Albert Gallatin, Secretary the Treasury, might well 
have been mentioned the section such projects, also the second ap- 
pendix. That report, which reviewed the status all existing public roads and 


canals under construction (or completed) the date was written, 
available. 


The Historical Register the United States, Part for 1814, edited Pal- 
mer, and published Palmer, Philadelphia, 
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DISCUSSION 


And finally, there seems reason believe that William Hammond Hall, 
California first state engineer, might have deserved little more attention. 
the West,” but possibly with too little explanation regarding his fine engineer- 
ing accomplishments. Irrigation received its first recognition state-wide 
basis the state California, and Mr. Hall was the first man have had 


FIG, 62.—WILLIAM HAMMOND HALL 


charge that project. skillfully did undertake such work that Captain 
Dutton “In this State (California), irrigationhas made greater 
advances than any other.” 

pioneer this field, was necessary for Mr. Hall develop stream- 
gaging methods and equipment well irrigation procedures. doubt 
was the nation’s foremost authority those subjects during his period ad- 
ministration. example his design work stream-gaging equipment 
the design the first American cable-way for making streamflow measure- 
ments. also improved the design two current meters, the first which 
his assistant, Grunsky described type.” The second current 
meter which made improvements was the Small-Haskell type. 
loaned several this second group the Geological Survey later 
years for conducting its work Colorado and Nevada. Surely, William Ham- 
mond Hall deserves having his picture this “Hall Fame” much his 
assistant, Grunsky. 


STEVENS, the writer’s opinion, this paper the hy- 
draulics classic the 20th century. not treatise hydrometry and 


First Annual Report the Irrigation Survey, (later converted into the 
Geological Survey). 


Cons. Engr. and Partner, Leupold and Stevens Instrument Co., Portland, Ore. 
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hydraulics. atreatise menwho sciences click the United 
States since the beginning its history. And this not all. Mr. Kolupaila 
preparing bibliography Hydrometry from the beginning writ- 
ten history, some 6,000 years ago, the present. Only one with his broad 
knowledge hydraulics and hydrometry and his versatility languages could 
successfully undertake such momentous task. 

example the author’s thoroughness, the writer aware the fol- 
lowing incident. has stated: “Some the sub-normal Nile 
flood stages were increased arbitrarily order unscrupulous rulers.” 
Mr. Kolupaila’s bibliography this nefarious desecration pinned down thus: 
“The rulers Egypt established land taxes called Kharag, corresponding 
the top the Nile flood that year. was the best interests the rulers 
have the Wafa (When the Nile reached the 16-cubit (28-foot) stage the sheikh, 
guardian the nilometer after the Arabs conquered Egypt, had proclaim 
the “Wafa,” beginning the annual fast.) proclaimed earlier and the top 
higher order entitle them start the tax collection earlier and 
increase the assessment. The scale was doctored for that purpose order 
the Caliph Omar: the cubits between and marks weré 
one-half. This fraudulent alteration resulted riots the population and 
the value the longest record river observation extending over 
1,300 years.” 

The writer wondering whether the relation forests water supplies 
could properly included the category hydrometry. the turn this 
century the forest-water supply relationship was one the major topics for 
the lecture platform and periodicals. The reason for this lies the fact that 
this was the time when the nation was debating, Congress and the press, the 
proposal that the government should acquire large areas forest lands and 
administer them National Forest Reserves. The conservationists were ur- 
ging their acquisition, but there was one hurdle congress had make before 
such course could legalized. 

Congress had authority over river navigation, but unless could shown 
that forests were substantial aid navigation congress was legally power- 
less acquire and administer forest reserves. The ardent conservationists 


therefore propagandized the nation with all the weapons the command. 
Two examples will cited: 


Gifford Pinchot, probably the greatest exponent forestry the world 
ever knew states:81 “The connection between forests and rivers like that 
between father and son. forests, rivers.” From this would appear 
that forestry substantial adjunct hydrometry. 

the crescendo direful prophecies, this one from the scrap book 
September, 1908, wins the blue ribbon: “When our forests are gone, the 


streams will dry up, the rivers will cease run, rain will fall more and 
America will desert.” 


One might add “and hydrometry will but memory.” One can imagine 
Teddy Roosevelt saying his tennis cabinet: “Go it, boys, the end justifies 
the means. Pull out all the stops. Make America forest conscious. Never 
mind the facts. We’ve got create forest reserves.” 


State Records for the Nile River,” Jarvis, Transactions, ASCE, 
Vol, 101, 1936, 1026. 


“The Fight for Conservation,” 1909, 53. 
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DISCUSSION 


One voices raised against falsification was lengthy 
and scholarly paper82 General Chittenden the U.S. Army Engineers. 

The writer has all known long time records stream flow which 
are analyzed with reference the forest cover their drainage areas. 

Thus, out welter conjecture and falsification, grew one our most 
valuable national assets. National forests now cover one-third our land 
area. With the sustained yield system now force not alone federal but 
also many private owners, are assured perpetual source lumber and 
corollary unrivalled recreation facilities, while soil erosion kept 
minimum those areas. The end did justify the means. 

One tie the forests have with hydrometry the fact that some forested 
areas are actually being cut off increase the water supply with the knowledge 
and acquiescence forest management. forestis great dissipater water 
evaporation from limbs and leaves and the extraction water from the 
soil growth processes. obvious, therefore, the forest removed 
more water will available for storage and use. appears, therefore, that 
forestry its many phases may definitely akin hydrometry. 


Transactions, ASCE, Vol. LXII, March, 1909, “Forests and Reservoirs Relation 
Stream Flow with Particular Reference Navigable Rivers.” 

Journal the Assoc. Engrg. Socs., Boston, July, 1913,“Forests and Their Ef- 
fect Climate, Water Supply and Soil.” 
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GENERALIZED DISTRIBUTION NETWORK HEAD LOSS 


Discussion Paul Constant, Jr., Marcel Bitoun, Claud Lomax, 
Joseph Maier and Thomas Miller, and Robertson 


PAUL CONSTANT, JR.4—The author describes mathematical method 
which pressure conditions existing water distribution system can 
obtained easily for different ratios Qp/ Before the method can used, 
there must available least one set measured data for two different 
ratios. Also, the various head losses are calculated under the assump- 
tion that they are result ratios where the loads are proportionately 
changed and all other system parameters remain unchanged. 

Although the restrictions imposed seem quite severe, the author’s method 
for calculating different head loss conditions useful for municipality distri- 
bution operations for “closed” system, that is, the operating end the busi- 
ness. The author very clearly indicates the usefulness his method for ob- 
taining intermediate pressure conditions. However, does not stress that 
application his method calculating head losses the operations 
end distribution problems. example where the author’s method may 
used effectively the determination pressure gradients system for 
“drought” year. 

Random checks the author’s calculations for and were made using 
the data from Table Anerror, An, greater than 2.46) was 
found. Further calculations were made determine the error, 
0.0344) due Aninn. This error was greater than 2%. These calcula- 
tions indicate that the author’s method within the limits engineering ac- 
curacy. fact, the errors are trivial because (a) estimates are made de- 
mands, (b) loads are grouped, and (c) error meter readings. 


The approximate value the relative error may determined 


easily using the general formula for errors, 


where 


which represents the relative error; the function several inde- 


Engrg. Dept., Midwest Research Inst., Kansas City, Mi. 
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represents errors the parameters. The general formula, Eq. can ap- 
plied the author’s equation, 


Qp | 
give the approximate relative error the calculation the function 


1.85 


which and are treated constants, and and are the errors 
the calculated values and respectively. 


The approximate relative error was found insignificant 


using and selecting and bethe errors previously calculated. 

There error the paper under the heading PROPORTIONAL LOAD 
CHARACTERISTICS WITH EQUALIZING STORAGE where the author solves 
for The number 0.67 should read 1.67. 

The author could enlarge upon the analysis ofa particular distribution sys- 
tem with the investigation flows the various pipelines. This natural 
consequence having the head losses. The flows can calculated simply 
means the Hazen-Williams hydraulic formula (constant density) the form, 


which the flow rate; the head-loss coefficient the pipeline; 
denotes the arbitrary constant account for flow units; and represents the 
head loss. 

The method calculations described the paper are little value the 
person concerned with improvements ofa particular distribution system; such 
as, (a) determination the best locationfor pipelines, (b) comparing alternate 
arrangements systems planned for construction, (c) optimizing system 
minimize the number sources required, (d) selection pipeline di- 
ameters for the best combination economy and performance. These prob- 
lems, along with others, seem the majority for consultants the field 
water and gas pipeline distribution analyses. This, part, was implied 
the author’s statement, “It must realized that these values and would 
changed either the load the piping the network Fig. were modi- 
fied any way, was changed.” 


MARCEL ASCE.—Publication, the author, the results 
the experiment analysis large, actual networks very valuable. Such data 


Chf., Design Branch, Div. Flood Control, Pennsylvania Dept. Forests and 
Waters, Harrisburg, Pa. 
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DISCUSSION 


are, unfortunately, too scarce. The author also complimented for hav- 
ing put evidence definite characteristics distribution systems with mixed 
types loading (residential and industrial). 

lengthy part the paper has been devoted the justification hy- 
draulic property proportionally loaded networks. Three distribution sys- 
tems increasing complexity were analyzed and used empirically justify 
the fact that “with proportional loads the percentage distribution flow in- 
dividual pipes balanced network constant irrespective the magnitude 
the total demand.” The author expressed this “contention” his con- 

This property has been used the writer and, presumed other hy- 
draulic engineers for the analysis balanced networks which the primary 
interest was the distribution discharges rather than the determination 
head losses. Actual flow values were replaced proportional figures, the 
largest being assumed 100 arbitrary units. Sometimes cumbersome fig- 
ures were thus eliminated. The writer wishes submit here demonstration 
the validity the hydraulic property empirically deduced Mr. McPherson. 

Consider distribution system without equalization storage governed 
boundary conditions expressed given ingoing and outgoing discharges. 
this system, let one balanced regime flow, characterized set 
ingoing and outgoing discharges, and the corresponding values flows the 
individual pipes that compose the system. 

have, for any junction this system, 


i=l 


which the number pipes connected the junction. 
any loop the system, have 


which the number pipes the loop. 


Let multiply all boundary discharges factor 0), and let 
the new regime flow the system. 


any loop Eq. can multiplied and yield 


+ 
} 
4 
| 
= 
th 
i=1 
4 
i=1 
are 
then 
i=1 
1 


a) 


Eq. multiplied yields 
i=1 
then 
i=1 


Eqs. and 12, valid for any junction and any loop respectively, show that; 


the regime hydraulically balanced, 


the flowfor every individual pipe equal the flowfor the same 
pipe multiplied the factor 


Moreover, this flow distribution the only one compatible with the given 
boundary discharges. Demonstration the unicity the flow distribution 
system where boundary discharges are given was furnished Angles 

Therefore, since the regime which all discharges are those charac- 
terizing multiplied bya constant factor, solution, and since this solution 
the only flow distribution the network corresponding the proportional 
load, have shown that the flow each individual pipe constant fraction 
the total load, which confirms Mr. McPherson’s findings. 

The number computations made the analyzer for the construction 
Fig. could therefore have been considerably reduced. single analysis for 
each network would have been sufficient provide one basic flow distribution. 


The second comment which offered here refers andn. For Mr. 
McPherson’s Eq. 


was stated that for =1, =K. 


However, two field tests were conducted for each location. fact, only one 
test necessary. The parameter actually the hydraulic resistance the 
system, and can computed from the physical characteristics the pipes 
which compose it, provided, course, that these characteristics (C, are 
known. Single field test will furnish the second equation necessary for the 
determination 

The last comment refers tothe relative influence non-proportional vari- 
ations the industrial loading and the residential loading the same net- 
work. felt the writer, although this not certainty, that deviations 


the unicity the solution network problems), Houille Blanche, No. May-June, 
1947, Grenoble, France, 209. 
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DISCUSSION 


from proportional loading within certain range will not result important 
differences the distribution flow the system. This would due 
“dampening” the difference throughout the system. Comments the author 
based onhis experience onthis point and especially results computations 
made the analyzer might prove conclusive. 


CLAUD LOMAX, ASCE.—The author fortunate having had 
opportunity compare actual flow conditions pipe network with the results 
obtained McIlroy Analyzer. The writer has not had that opportunity and 
must present his comment the light this lack experience. 

From theoretical standpoint, the technique proportional load hydraulics 
can firmly supported. The comparisons that the author makes between the 
Analyzer results and the actual flow distributions and the pressures, 
show some discrepancies. The differences may attributed, part, the 
following: metering limitations, Fluistor selection limitations, and mistakes 
that may have occurred the Analyzer solution. The conformity 
the comparisons effected the reliability the field measurements flow 
distribution, pressure losses, and data such pipe sizes, lengths, and friction 
factors. The writer the opinion that the results presented this paper 
are well within the range discrepancy expected. 

the light the writer’s experience, which has been primarily with dis- 
tribution networks for small western cities, word caution general ap- 
plication proportional load hydraulics order. the smaller cities the 
preparation skeletonized arterial network usually not feasible, be- 
cause there may onlyone run large pipe from pumping station res- 
ervoir, and the remainder the network made in., in., and 8-in. 
sizes. The author stated, “This the smallest distribution district Phila- 
delphia and the only one which fire flows governed the design future ar- 
terial piping requirements.” For smaller cities the fire flows are usually 
controlling. For western areas irrigation loads are significant parts the 
total system load. Consequently, the writer’s opinion that application 
proportional loads for design conditions flow networks small cities not 
feasible. 

The author has not stated, but implied that the source magnitudes 
must proportional wellas the load magnitudes. The smaller cities draw 
from both reservoir and well sources, and usually all sources are required 
satisfy the maximum demand condition, whereas conditions less than maximum 
require reservoirs that scattered throughout the community. 

would very difficult achieve proportionality source magnitudes 
system where reservoir water surface elevations and pump discharges are 
fixed. One would need additional pump capacities each pumping station, 
which were proportional the percentage contributed that station well 
proportional the total load increase. Fire loads, when superimposed 
the maximum day demand for irrigation, domestic, and industrial usage, 
not increase the load the system proportionately since these fire loads are 
placed specific hydrants selected throughout the community. The use 
proportional load hydraulics the smaller cities may valuable, judiciously 
applied, determine the capacity the system for flow distributions greater 
than the maximum day demands. clue the behavior the system with 


Assoc. Hydr. Engr., Albrook Hydr. Lab., Washington State Univ., Pullman, 
Wash, 
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fire loads can obtained this way. However, the pressures and flow dis- 
tributions inthe immediate vicinity those hydrants being used may greatly 
error. 

Having used the Analyzer for considerable number flow dis- 
tributions, the writer, and probably others having had close contact with sev- 
eral Analyzer distributions, impressed the feeling one gets for 
“healthy” network distribution system. well-designed system the omis- 
sion section pipe, other than the critical section leading reservoir 
pumping station, usually does limit the capacity the system. 
Through errors wiring selection incorrect Fluistors, one accidentally 
sees this demonstrated network, that not really obvious that 
error has been made until check procedure discovers this error. other 
networks, which are “weak,” the omission almost any pipe the network 
distribution system will reduce its capacity. Proportional load hydraulics 
probably much more applicable “healthy” system than “weak” one. 

With very low flows, conditions present wherein laminar flow exists 
some the pipes. apply proportional load hydraulics and increase the 
flow the point where these flows become turbulent, certainly unwarranted. 

Referring the author’s gas flows, should noted that gas 
flow distributions are related the 1.85 power the flow when the Fritzche 
Panhandle formulas are used. 

The author tobe complimented for having presented, with supporting evi- 
dence its validity, technique analysis which, judiciously applied, will 
serve designers flow distribution systems. should obvious that one 
does not require Analyzer analysis nor complete system analysis 
initiate use this technique. Having pressures critical points dis- 
tribution network andassuming the system proportionally loaded, one can then 


evaluate the total capacity the system. course, one can start with 
network analysis well. 


JOSEPH and THOMAS MILLER,® A.M. ASCE.—The propor- 
tional loading principle, presented this paper, does not seem ap- 
plicable gas system networks. The term “proportional load” described 
the typical design assumption that each individual consolidated load varies 
fluctuates about its mean direct proportion the total system fluctua- 
tions. Thus the maximum and minimum hour demands the average and 
maximum day would taken fixed multiples, not only the annual aver- 
ages for the district whole, but also for each individual local load. has 
been the writers’ experience that consolidated load fluctuates according 
the characteristics the primary gas load which composed. From 
this load characteristic fluctuation, can further concluded that fixed 
multiply representing total system demand variation apparently not appli- 
cable individual consolidated loads gas distribution system. 

The composition ofa consolidated gas load will include gas used for one 
more the following purposes: commercial, industrial, house heating, cook- 
ing, water heating, refrigeration, and drying. consolidated load consisting 
industrial customers will not realize appreciable yearly variation, while 
consolidated load composed residential house heating customers will fully 


Asst. Engr. Distribution Design, Philadelphia Gas Works, Philadelphia, Pa. 
Staff Engr., Philadelphia Gas Works, Philadelphia, Pa. 
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DISCUSSION 


affected temperature and wind velocity variations. can seen that each 
the above mentioned gas usages has specific characteristics which when 
considered network analysis cannot represented fixed multiply. 
addition, seasonal industrial customers are supplied during the non-heating 
season. 

illustrate gas load fluctuation, letus use and Fig. sche- 
matic flow breakdown winter day, which the total hourly sendout, 
refers the industrial customers, indicates the residential non-heating 
customers, and the residential heating customers, while Fig. which 


WINTER 


SUMMER 


FIG, 


the residential non-heat customers plus seasonal industrial customers, 
and refers the non heat (base load) residential heating customers, 
represents the same network summer day. Although the summer hour 
sendout 60% the maximum winter hour sendout, the consolidated loads 


the summer day Qo, are not 60% fixed multiplies the corresponding 
maximum winter loads. 


ROBERTSON, seems desirable point out that the 
solution pipe network flow and head loss problems with digital computer 
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after the Hardy Cross technique has been improperly termed relaxation. Al- 
though many authors use the term relaxation include the iteration pro- 
cess employed the Cross approach, there distinction that should ob- 
served. The error terminology not surprising since the introduction 
the iteration technique obtaining solutions such problems 
was but his work where did introduce 
the technique (as noted relaxation technique was proba- 
bly used early 1823 Gauss for solving the simultaneous equa- 
tions encountered surveying) later named relaxation Southwell. 

The pipe network analysis technique sometimes termed the Cross-Doland 
method, since was the suggestion that Cross considered 
the application the methods numerical successive approximation net- 
work problems. noted recently pipe networks may 
solved either iteration relaxation methods. The relaxation method 
primarily limited paper and pencil type calculations, which ideally 
suited and for which was originally devised, clearly stated 

spite the fact that relaxation techniques are applicable various types 
fluid flow problems, which may now attempt solve with the helpof high 
speed digital computers, solutions such problems digital computers are 
rarely, ever, relaxation solutions. They are merely iteration solutions—the 
two terms are not synonymous. This because relaxation employs judgment 
and reasoning, not intuition, every step. These are beyond the usual men- 
tality whichis programmed into digital computer. Sucha machine can merely 
plug along, very high rate speed course, making one simple calcu- 
lation after another accord with rigid prescribed program. Relaxation 
techniques, the other hand, accelerate the iterative convergence solutions 
through various Thus, computation first applied regions 
where the divergence proper local result greatest. addition, such 
techniques block and over relaxation are employed where suitable. Since 
judgment required deciding when and what technique employ this can 
hardly programmed into dumb computer. Or, can be, the program 
just too complicated, too involved. Quite probably nothing impossible ona 
computer, but better simple computation many times than try 
complicated one that can achieve the same computation one cycle. Thus, 
Southwell himself has said, best set computing machine work with 
iterative process. the present writer’s understanding that such in- 
deed the case with distribution network analyses. 

The question the value the exponent the power-law (sometimes 
incorrectly termed exponential law) head-loss relation raised author. 
Although this exponent may assume values between 1.75 and 2.00, network 


Illinois, Engrg. Experiment Sta. Bulletin No. 286, November 13, 1936. 

“Gauss Gerling Relaxation,” Forsythe, (translation and commen- 
tary letter Gauss dated December 25, 1823) Mathematical Tables and Other Aids 
Computation, Natl. Research Council, Vol. 1951, pp. 255-258, 

“Simplified Analysis Flow Water Distribution Systems,” En- 
gineering News-Record, Vol. 117, October 10, 1936, pp. 475-477. 

“Relaxation Methods for Pipe Networks,” Chapman, Civil Engineering 
and Public Works Review, Vol. November, 1956, pp. 1247-1248. 

“What are Relaxation Methods,” Forsythe, Chapter Modern Math- 
ematics for Engineers, Beckenbach, McGraw-Hill Book Co., 1956. 
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analysis progresses best single value used for any given analysis. 
Certainly the author’s proportional load hydraulics method not applicable 
several pipes different values occur the system. Network solutions 
techniques may not limited fixed-m (versus variable several sys- 
tems, but the exigencies the approaches the solution, such are being 
discussed, often limit them. Certainly any paper and pencil analysis (plus 
slide rule desk calculator and table powers numbers) whether merely 
iterative relaxation, can easily employ different values for each 
The network analyzer may used with other than 1.85. Different 
values could even employed for different pipes, but this rather trouble- 
some. Digital computer iteration solutions have been employed for values 
2.00 and less, but with the same value for all pipes any given network. 
Henry and note, however, that the use non-integer val- 
ues the exponent tied large partof IBM 650 machine for the process 
raising the flow rate this power. cannot this nearly 
one can with analog computer (such slide rule, for in- 
stance). This extra effort reduces the number pipes that can handed over 
that number which can handed when 

The programming complexities using different values for the several 
pipes system are obviously much greater, since special subroutine 
would needed for each value employed. Either extensive table 
such results for limited set values would needed some analog 
auxilliary must developed which the digital computer may refer quickly 
when needed. The first method probably impractical due the ac- 
cess time required, while the second not available. Since the range flows 
given pipe usually small, there seems need consider pipe 
lines variable would thus appear that for some time come our so- 
lutions network systems limited few fixed values the head loss ex- 
ponent. Finally, course, there remains the practical question whether the 
exponent can determined accurately enough justify considering anything 
more than average value. 


“Computer Evaluates Design Gas Gathering System,” Henry and 
Peritz, Oil and Gas Journal, Vol. 57, No. 17, April, 20, 1959. Also “Evaluation 
Natural-Gas System Design with Electronic Computer,” ASME Paper 58-Pet-28 
(for condensed version see Mechanical Engineering, June 1959, pp. 58-59.) 
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THE FOURTH ROOT n-f 
Discussion Gerald Lacey and Rajaratnam 


GERALD author writing his paper has, manifestly, two 
objects mind: first the presentation new friction-factor diagram for 
the use practical engineers, and secondly that challenging the validity 
logarithmic flow formulas. 

Provided that the engineer makes all the information his dis- 
posal would seem that makes little difference whether uses the Moody 
Diagram which effectively that Colebrook and White, the Blench dia- 
gram. The author suggests that the latter should employed and judged 
results. There little doubt that the author’s diagram can usefully em- 
ployed but such use can neither prove not disprove the validity the loga- 
rithmic flow formulas which, present, hold the field. 

Those who have made recent publications the theory flow 
pipes might well inclined conclude that, far closed conduits are 
concerned, theory now approaching finality. The use Colebrook the 
title “Smooth and Rough Pipe Laws” supports this view. these laws, so- 
called, their logarithmic form which the author believes lack completely 
logical foundation verification the range practical data. Such “proofs,” 
more properly physical explanations, have been put forward the aca- 
demic world, regards patently illogical and mathematically obscured. 
The author puts forward alternative his Fourth Root Theory based the 
postulate universal flow formula for all boundary types which, would 
appear, would embrace open channels with moving boundaries well pipes 
with rigid boundaries. The author’s diagram, practical although may be, 
essentially vehicle for the presentation the author’s theories flow. They 
not lack originality. 

The author draws attention the fact that Prandtl was careful explain 
that his proof logarithmic velocity distribution did not apply pipes but 
would have been value the author had stated when, ever, the proof 
was valid—in the case open channel for example. The point relevant 
has hesitation making use observations open irrigation channels, 
transporting sediment, and with moving boundaries. 

rigid boundaries are concerned that with full 
turbulence and low values the relative roughness the logarithmic equation 
can replaced over large range the Manning formula with little error. 


Cons., Sir MacDonald and Partners, London, England. 
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The contention the author that “believed” that Manning’s slightly 
dependent the size pipe debatable. 

When the author refers “absolute roughness,” and the well known fact 
irrigation canal practice that Manning’s for exact design has reduced 
for large channels the same apparent boundary roughness, stands 


ground. The writer more than ago (about 1930) was responsible for the 
formula which the author relies! 


(this equation being metric units), and commented that 


“the coefficient was one absolute rugosity, equal that Kutter and 
when the hydraulic mean depth was metre, but remaining 
true constant long the size the silt particle was constant. 


The need vary the rugosity coefficient with the size the channel 
vanished.” 


This equation however applies open channels with moving boundaries, and 
possibly very different velocity distribution from that channel witha 
rigid boundary. 

There doubt that the author’s diagram has the merits great simpli- 
city. Whether the practical engineer should weaned from his predilection 
for the Manning equation when dealing with pipe flow and should use the fourth 
root equation instead openquestion. The use the fourth root throughout 
may appeal tidy mind. 

When the author states that content treat his value code 
number intuitively observer recognizes faces,” will find 
many who will agree with him. The difficulty measuring roughness quanti- 
tatively stumbling block and this matter judgment and experience are 
value. This may some comfort the practical engineer but would 


like some assurance towhether should use the author’s “e” the more 
usual “k.” 


new method solution for the problems flow resistance pipes. 

The writer wants make few observations the projected universality 
the fourth root n-f diagram. Firstly, the author seems believe that the 
Blasius law, given the equation 


describes the turbulent flow resistance smooth pipes even for values the 
Reynolds number higher than 10°. But widely known that for values the 
Reynolds number higher than the flow resistance better approximated 
the Prandtl-Karman equation 


Proceedings, Inst., Civ. Engrs., Vol. 229, pp. 372, 373. 
“Research Fellow, Civ. and Hydr., Indian Inst. Science, Bangalore, 
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fact, even the lower range, Eq. agrees with Eq. quite 
that when single law resistance desired, Eq. that normally 
adopted. 

Secondly, regarding the universal adaptability the fourth root n-f equa- 
tion, the author has suggested that if, the equation 


1/4 


the thickness the laminar sublayer (which believed exist near the 
smooth boundary for turbulent flow) substituted for the Blasius formula 
will obtained. The thickness the laminar sublayer the 
equation 


which the thickness the laminar sublayer, the shear velocity, 
and the kinematic viscosity the fluid. 


Substituting 


Substituting Eq. Eq. 


against constant which the Blasius law. 
The derivation!§ the equation 
1/4 


equating the values for Mannings, Nikuradse and Fourth Root formulas 


150 not clear. Further, the modification Eq. into the equation 


“Boundary Layer Theory,” Schlichting, pp. 402, 413. 

“Unification Flow Formulas,” Blench, Proceedings, Vol. 
Lisbon, 1957, 


an 
2 
= 
which the mean velocity flow the Eq. 11. 
V8 
2 


July, 1960 


seems rob any possible physical significance. Eq. the outset 
looks like one the many exponential formulae proposed different times. 
For example, recent times, has proposed formula that can 


shown 
1/6 


which (Colebrook’s) equivalent sand roughness. 
Based Eq. 18, can shown that 


which agrees closely with the formula proposed the writer for open channel 
flow. 

generally known that decreases with increase the size the pipe 
recently confirmed the careful experiments Webster and 
how the author made independent the size the pipe the multiplication 
factor Webster and Metcalf data does not confirm this shown below. 


0.0246 3.0 0.02695 
0.0244 5.0 0.02790 
0.0236 7.0 0.02770 


The equation for the transition line rough turbulent flow has been given 


R1/4 0.80 


has been that the Nikuradse’s sand roughness and the 


thickness the laminar sublayer, the transition starts <4.0 and ends 
6.0. Based finding, the equation the transition line the Moody 


Chart has been shown 


This equation has been shown hold good many cases excepting few dis- 
crepancies recently reported. Since, this present case, seems give 


“Recent Developments the Flow through Pipes and New Method Approach 
for Rough Turbulent Flow through Pipes,” Advani, Irrigation and Power, Octo- 
ber, 1957. 

Contribution Turbulent Flow Open Channels,” Govinda Rao and 
Rajaratnam, Presented the 30th Research Committee Meeting the P., 
India, June, 1960. 

“Friction Factor Corrugated Metal Pipe,” Webster and Met- 
calf, Proceedings, ASCE, Vol. 85, No. September, 1959. 
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direct magnitude the roughness projections, the same line argument may 
not valid and hence the validity Eq. doubtful. 

agreed that the Colebrook-White formula, which the curves the 
transition range are based, does not give dependable for all surfaces. 
fact Colebrook29 himself said that his equation represents mean curve for 
surfaces “tar-coated cast iron, wrought iron and galvanised iron.” has 
been claimed Morris2! that the concept the “equivalent sand roughness” 
not sound because Nikuradse’s data based gives rising curve 
the transition region against falling curve predicted the Colebrook- 
White equation. Further, the Colebrook-White equation does not make any 
provision for conduit surfaces having either rising horizontal charac- 
teristic the transition region. Thus clear that Moody’s Chart does not 
give dependable values for all surfaces the transition region. But this 
difficulty still there the fourth root n-f diagram also. Hence new ra- 
tional method approach should sought, having eye for the physical pic- 
ture the flow. this direction, admirable attempt has recently been 

The writer suggests that cases where the flow the wake-interference- 
type the case many surfaces open flow, may advisable 
retain the concept equivalent sand roughness until better solution 
evolved. 

Regarding problem which was worked out the author, normally 
that due the additional effects the free surface, the existence 
secondary currents and non-uniform boundary shear, the equations resis- 
tance for open channel flow are likely different from those for pipe flow. 
Based the theoretical equations Garbis Keulegan and the experimental 
results Mr. Powell, the has formulated set rational formulas 
for flow resistance open channels with new conversion formula for and 
and another formula for the transition line rough turbulent flow. 

Acknowledgment.—The writer grateful Govinda Rao for his 
interest the preparation this discussion. 


“Turbulent Flow Pipes with Particular Reference the Transition Region Be- 


tween the Smooth and Rough Pipe Laws,” Colebrook, Proceedings, Inst., 
Civ. Engrs., London, 1933, pp. 133-156. 


«Design Methods for Flow Rough Conduits,” Morris,” Proceedings, 
ASCE, Vol. 85, No. July, 1959. 


“Flow Rough Conduits,” Morris, Transactions, ASCE, Vol. 120, 1955. 
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CONSERVANCY DISTRICTS FLOOD CONTROL 


Discussion Henry Tebow 


HENRY ASCE.—The author has presented refreshing idea 
for developing additional local participation and partnership solution 
problems “in which Federal interest predominant,” such the problems 
major rivers. 

The conservancy district approach invaluable tool the consultant 
who sometimes has difficulty, due special interests and their desires, 
winning acceptance over-all plan, which will best fit the district 
whole. 

There should rewarding work from conservancy district operation, its 
consultants, and the state and federal engineers charged with coordination 
planning, design, and construction. The conservancy district will make prac- 
ticable the application “New Approach Local Flood sorely 
needed economize flood control requirements for the nation. 

The Colorado Conservancy District Act May, 1937, specifically paved the 


way for the construction the Colorado-Big Thompson Project, Federal 
water project designed divert surplus waters the western side the 
Continental Divide the “water-short” acres the eastern slope Colorado. 
Though repayment planning for this Federal project allowed credit for flood 
control, operation the reservoirs maximize the storage snow-melt 


floods does provide flood control. The history events for the project 
follows: 


Initial Settlement and Irrigation Area 1859 
First State Survey diversion ($25,000) 1889 
First Bureau Reclamation investigation 1904 
Reservation for Reclamation use Rocky 

Mountain National Park Act Congress 1915 
Consulting Engineer Report 1933 
Public Works Administration Funds 1935 

Conservancy District Organized 1937 
Project authorized (Senate Document 80) 1937 
Initial construction (replacement dam) 1938 
First power produced (To Public Service Company 

Colorado) 1943 
First supplemental water delivered 1947 
Initial year major supply lands 1954 


April, 1960, Cloyde Chambers. 
Genl. Engr., Design and Constr. Div., Region Bur. Reclamation, Denver, Colo. 


“New Approach Local Flood Problems,” Herbert Vogel, Proceedings, 
ASCE, Vol. 86, No. January, 1960, pp. 
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Completion for Supplemental Irrigation 1956 
Project Completion (Final Powerplant six) 1959 


Project Data 
Supplemental irrigation for 720,000 acres. 
Power developed for pumping and sale 183,950 kw. 
Average surplus water diverted 257,000 Acre Feet. 
Average annual power from the water 3/4 billion KWH. 
Annual income Treasury approximately million. 
Cost $159 million—all repaid. 
Power costs repaid with interest. 
Project water value 1954 $22 $41 million crop. 


The Colorado Conservancy District Act 1937 provides for the organiza- 
tion conservancy districts any District Court upon petition required 
number property owners. The court appoints board directors that has 
the power appoint officers, acquire and hold property, appropriate 
water, enter into contracts, levy taxes and assessments, allot water, 
and generally administer the affairs the District.4 

reflection the indirect benefits accruling tothe entire area included 
conservancy district, the act authorizes assessment all real and 
personal property not exceed one-half mill during the construction period 
nor exceed one mill thereafter, except case default deficiency, when 
additional one-half mill permitted.4 

The act further provides procedures for allotting water individuals, 
municipalities, and irrigation districts and the collecting special taxes for 
the use such water. The board has the power contract with the United 
States government for construction, when authorized election property 
owners. 

The historical record events for the Colorado-Big Thompson Project 
shows that the Conservancy District Act 1937 was effective initiating 
plan construction just about after initial surveys were started. 

Western pioneers, starting with the Union Colony and early irrigation near 
Greeley, Colorado (named after Horace Greeley—author the quote, “Go west 
young man”) overestimated the river water supply. Late summer and every 
year drought periled the very settlements the western end the rail- 
roads. The more without construction water plan were certainly 
due lack conservancy district law. 

Native American conservatism, and perhaps difficult basic communication, 
can given the reason that lapsed before the Miami Conservancy 
District idea went west from Ohio. 

Not with regard the hydraulics the MCD dry dams. Ivan Houck 
and the fine engineering papers hydraulics, published MCD, came west 
the United States Bureau Reclamation the 1930’s. The legal values 
Miami Conservancy District were slower moving westward. 


reported the Tech. Record Design and Const., Vol. Colorado-Big 
Thompson Project, Bur. Reclamation, Dept. the Interior, 1957. 
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118 and 119. These pages should placed immediately after 123 
and before 124. 


123. Please insert the words 
radius 0.056 was omitted. With W/D 1.5 and offset 0.05 
Fig. 12, the pressure parameter for point downstream 
just prior the words 
from the P.I. 
line the last paragraph. 
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